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Chaos may make black holes bright 
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Black holes cannot be seen directly since they absorb light and emit none, the very quality which 
earned them their name. We suggest that black holes may be seen indirectly through a chaotic 
defocusing of light. A black hole can capture light from a luminous companion in chaotic orbits 
before scattering the light in random directions. To a distant observer, the black hole would appear 
to light up. If the companion were a bright radio pulsar, this estimate suggests the black hole echo 
could be detectible. 
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Black holes evade direct detection precisely because 
they are black. The existence of black holes hidden be- 
hind accretion disks or in the centers of galaxies have 
been inferred from astrophysical observations. Despite 
these indirect observations, we cannot know for certain 
that the compact objects lurking there are in fact Ein- 
stein's black holes. Any detection which can see in very 
near to the event horizon would provide more incriminat- 
ing evidence for their existence. In this Letter we describe 
how chaotic scattering of light in an inner regime around 
the event horizon could effectively render the black hole 
bright. If perturbed, a stochastic region develops around 
the last unstable photon orbit. The disturbance could be 
an orbiting companion or the emission of gravitational 
waves or any asymmetry in the evolution. The black hole 
can then trap incident photons in the stochastic region 
for some time before throwing off half and absorbing the 
other half, effectively shrouding it in light. We estimate 
the cross-section for this gravitational defocusing gener- 
ically and then illustrate with an extremal binary black 
hole spacetime. 

Around an isolated black hole of mass M and charge 
Q, light follows the simple orbits 

The motion of massless particles depends only on the 
impact parameter b = L / E and not on the energy E and 
the angular momentum L separately. Given the mass 
and charge of the black hole, there is a critical value of 
the impact parameter, b c , at which light gets trapped in 
a perfectly circular, unstable orbit. For smaller impact 
parameters, light will fall into the black hole while for 
larger impact parameters, light escapes. Just above b Cl 
the light can come close to the circular orbit executing 
one or more full rotations before being cast off. Only 
black holes are compact enough to bend light by more 
than 7r. The phenomenon of back scattering by a full 
7r is known as the glory 0] and was thought to be the 
weakest way to observe black holes. 

In the stochastic region around a black hole pair, the 
last unstable orbit becomes the site for chaotic scattering. 
The lone periodic orbit is replaced by a glut of periodic 
orbits. These proliferating orbits are packed so densely 
into phase space that they form a fractal set. Fractals 
are a way of maximizing the area while maintaining a 
bounded volume. Light scatters chaotically as it skips 
from one periodic orbit to another. The cross-section 
for catching light in multiple windings around the hole 
is then amplified. As well, the black hole hangs on to 
the light for longer and reemits the light more evenly. 
A bright light directed onto the black hole, say from a 
pulsar companion, could illuminate the black hole for a 
time before the light decayed away and the star fell dark 
again. 

The defocusing cross-section can be approximated gen- 
erally in terms of the topological features of the chaotic 




black hole 



0^=271^ A b 

FIG. 1. A black hole with a bright companion captures 
light around R — 3M for a few windings before the light 
scatters off. 

fractal set. The importance of the general approach is 
that the scattering cross-section for any candidate system 
can be estimated in this manner. We assume isotropy 
which will in fact be broken by an actual perturber. The 
perturbation is largest at the point of closest approach 
and is presumably larger in the orbital plane. This whole 
system is also moving relative to the observer. For the 
purposes of estimating the magnitude of defocusing, it is 
reasonable to assume isotropy. In fact, isotropy could be 
approached for complex orbital motions since light may 
no longer be confined to a plane. The cross-section a is 
roughly the geometric area of the annulus 

a = 2nb c Ab (1.2) 

as shown in fig. [j]. As light is shot at the black hole with 
impact parameter near 6 C , it will travel on nearly peri- 
odic orbits for a time before diverging from these unsta- 
ble worldlines. The cross-sectional thickness A6 is thus 
given by the number of periodic orbits times the thick- 
ness around each orbit in phase space e. With N(n) the 
number of fixed points lying on periodic orbits which ex- 
ecute n windings around the black hole, a can be written 
as a sum over all winding numbers, 

oc 

o = 2^blY^N{n)e{n) . (1.3) 

n=l 

The number of fixed points is given by the topological 
entropy 

S= lim -N(n) (1.4) 

n — >oo TL 

so that N(n) ~ e Sn in the limit of long orbits. Pertur- 
batively, the deviation is Sr/r ~ A with A the Lyapunov 
exponent, a measure of the instability of the orbit. The 
width in phase space is then e(n) ~ be~ Xn The instabil- 
ity can vary from orbit to orbit. As an approximation 
we take A to be the average over all the fixed points. 
Although the Lyapunov exponent is a notoriously coor- 
dinate dependent quantity, in this context we are measur- 
ing the instability in units of windings around the black 
hole and the winding number does not vary from observer 
to observer. Under the assumption of ergodicity, half of 
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the photons fall into the black hole and half are cast off, 
so we divide the cross-section in half to find 




= ^£(e X p(S-A))" 



71=1 



Performing the sum gives 



= irbi 



1 - e( 5 - A ) 



(1.5) 



(1.6) 



We can recast a in terms of the fractal properties of the 
set of periodic orbits. The fractal dimension is defined as 



D = lim 



lniV(e) 
^oln(l/e) 



(1.7) 



where N(e) is the number of boxes of size e needed to 
cover the set. We can relate D to the entropy by noting 
that e oc e~ An for orbits of length n, from which it follows 
that 



D = lim 

n— >oc 



lniV(e) 
nX 



(1.8) 



Again, this assumes A is the same across the set or that 
a suitable average will fare well enough. There is an 
entire spectrum of weighted entropies and dimensions to 
cha ract erized the fractal for which similar rela tion ship s 
to (1.8) have been conjectured Using (1 

we then estimate the cross-section to be 



in (1.6) 
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e -A(l-D) 

I _ e -A(i-D) 



(1.9) 



We could have deduced the cross-section directly from 
the fractal property. The proliferation of periodic orbits 
form a fractal set which fills an area in the oribital plane. 
The set, by definition of the fractal dimensionality, there- 
fore has a thickness. That thickness is determined by the 
fractal dimension. According to the usual determination 
of the width of a fractal set, the number of fixed points 
in the direction dr which can be covered by boxes of size 
e is given by N ~ (e/b)~ D . The length of this set is then 
Ne = b(e/b) 1 ^ D with e ~ be~ Xn . Sum min g and dividing 
by 2 we again derive the area of eqn. (1.9). 

We can conservatively evaluate a by using the unper- 
turbed values. The dimension of the boundary is zero 
in the unpertu rbed , nonchaotic system. To estimate A, 
we vary eqn. ([hi]) around the unstable circular orbit. 
For a chargeless Schwarzschild black hole of mass M, 
the last unstable circular orbit lies at R = 3Af with 
b c = 3V3M. The perturbed radial motion grows as 
SR/R ~ e™*, or in units of <f> = 2irn, the exponent is 
A = 2?r. We then estimate a - 0(lCT 3 )7r&;;. For an 
extremal black hole (Q — M), the last unstable circular 
orbit occurs at R = 2M with b c — AM. At second order 
we find SR/R ~ e"^/^. Measuring <j> in terms of the 
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FIG. 2. The color coded fractal basin boundaries for two 
extremal black holes with Q = M — 1. The initial location of 
the photon is painted black if it sticks to the mass at (0, 1), 
grey if it sticks to (0, —1), and white if it escapes. The initial 
radial velocity is r = and </> > is set by eqn. (1.15). The 
error tolerance was ~ 10 -8 . 
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winding number, <f) ~ 2irn, we read off A ~ yl-n. For 
an extremal black hole, the cross-section is larger with 
a ~ 0(lO _2 )7r6 2 ; that is, a hundreth the capture cross- 
section. 

Unlike other more conventional estimates of the width 
of a stochastic layer, this estimate requires knowledge 
of only a few simple properties of the system and does 
not require a complicated examination of the dynami- 
cal equations. Unlike other more conventional estimates, 
there are inherent shortcomings. As with the glory cal- 
culations |l]], the effect is dominated by the trajectories 
which wind around the black hole the fewest number of 
times. These trajectories may not model the fractal set 
of underlying periodic orbits as well as those which ex- 
ecute many windings and spend the most time chaoti- 
cally s catt ering off the set. In other words, the sum in 
eqn. (1.5) is dominated by the first few n while eqn. 
(1.4) is a large n limit. Given this caveat, we can see 
how well the approximation fares in a given dynamical 
system. To explicitly illustrate, we consider scattering 
around two extremal black holes. A pair of black holes 
with equal charge and mass are able to coexist in a static 
configuration with the electrostatic repulsion caused by 
their charge just balancing the gravitational attraction of 
their masses. While the resultant Majumdar-Papapetrou 
spacetime is static, the geodesic flows around the 
pair of black holes are known to be chaotic |||7| . We iso- 
late fractal basin boundaries |f7|-|i"l"| for massless particles 
as has already been done for massive particles J7j . 

The Lagrangian for motion in this space can be written 
in isotropic coordinates as 



while two are dynamical 



c = --u- 2 P + 
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with (l 2 — sin 2 99+ <j) 2 and an overdot denotes differenti- 
ation with respect to an affine parameter. Schwarzschild 
coordinates are recovered with R = Ur. The metric com- 
ponents are determined by 



rr , M M 

U = 1 + — + — 

ri r 2 



(1.11) 



where r\ is the coordinate distance from one black hole 
and r 2 is the coordinate distance from the other. The 
event horizons occur at n = and at r 2 = 0. Wc 
place a black hole with mass M and charge Q = M at 
(r,0,(f>) = (M, it/2, 7r/2) and an identical companion at 
(M, 7r/2,-7r/2) so that 



r 2 - 2Mr sin + M 2 
r 2 + 2Mrsm<b + M 2 



(1.12) 



The geodesic motion is found by evolving the conjugate 
momenta, Tl q = dC/dq. Two coordinates are automati- 
cally conserved with 
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n t = — = e 
n fl 
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= r 2 U 2 9 = Lg 



IU = r 2 U 2 ^ 
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U 2 r 
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and evolve according to the equations tl q = dC/dq. 
There is an additional constraint equation, equivalent to 
the conservation of energy obtained by setting C — so 
that the photons travel along null geodesies: 



f 2 + r 2 6 2 = E 2 



(1.15) 



We have assumed that Lg = and so consider motion in 
the plane defined by the binary system (9 = tt/2). 

We numerically evolve the goedesic equations, en surin g 
that the energy is conserved according to eqn. ( 1.15 ). 
In isotropic coordinates, the last unstable circular orbit 
around an isolated black hole occurs at r = M, So, 
in the black hole pair, we expect to see chaos around 
n = M and r 2 = M, which is precisely what we find 
as illustrated by the fractal basin boundaries of fig. || 
|j],H ■ We look at an initial slice through the plane of the 
binary system. We submerge the black holes and their 
surrounding area in a bath of light, a photon at each 
locati on in space with r = initially and (f> > set by 
eqn. ( |1.15| ). We color code the initial location black if the 
photon which originated there fell into the hole at r\ = 0, 
grey if it fell into r 2 — 0, and white if it escaped. The 
result is the mixed, fractal basin boundary as shown in 
fig. The lower panel focuses in on one region to show 
the repetition of the fractal structure. The dimension of 
this boundary is 1 + D ~ 1.1. 

If we take Ab to be roughly the size of the white strip in 
the lower panel of fig. [|, we would estimate Sr/M ~ 0.02. 
Of course, the meaning of this coordinate thickness is 
ambiguous. Still it is reassuring that it is comparable in 
magnitude to the value we would have guessed from our 
generic approximation using the unperturbed A and the 
measured value of D ~ 0.1, which gives Ab/b c ~ 0.02, 
although the extreme agreement is certainly fortuitous. 
The cross section is again about a ~ O(10 _2 )7r& 2 . 

It is not immediately obvious if the approximation un- 
derestimates or overestimates the corss-section. Some 
elements are underestimated. For instance, there are ad- 
ditional contributions to the defocused light from the in- 
ner regions of the binary sys tem which are not included 
in the estimate of eqn. (l.E). Also, any realistic astro- 
physical black hole with a companion will only be more 
chaotic than this crisp example. The metric will not be 
static and it is likely that the evolution of the spacetime 
itself will be chaotic. The companion needed to provide 
the luminous source may well be on an unpredictable 
trajectory contributing to the probability of defocusing. 
On the other hand, the calculation was restricted to the 
orbital plane where the effect would be largest and the 
relative orientation of the observer may decrease the sig- 
nal. 
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FIG. 3. A black hole, pulsar, Earth configuration. 

We can use the estimate of the cross-section from the 
Majumdar-Papapetrou spacetime to give rough predic- 
tions for the observability of the effect. To a distant ob- 
server, a black hole with a pulsar companion will appear 
to radiate with a luminosity 

£bh — i* ( — - — ( — ) (1.16) 

where L» is the observed luminosity of the star, the bi- 
nary lies at a distance d from the Earth, and d* is the 
relative star and black hole separation as shown in Fig. 0. 
The cross-section of the beam is roughly ndl sin 2 (Obeam) 
where Bbeam is the half-angle subtended by the pulse. 
We further assume that the time it takes the beam to 
sweep over the black hole r* is comparable to the char- 
acteristic time t for the captured light to decay off the 
black hole. We take the companion to be at a distance of 
- 10"A/ B h- Using a ~ 0(l(T 3 )7r27A/ 2 and assuming 
Obeam ~ 5° we estimate 

L BH > 3.6 x l(T 2n+1 L m . (1.17) 

Notice that this is the most pessimistic estimate. The 
black hole luminosity is suppressed relative to the inci- 
dent luminosity by a /abeam- Since abeam is much larger 
than the scattering cross-section, the returned radiation 
looks small. However, the incident radiation can accu- 
mulate as the pulse returns again and again feuling a 
larger returned radiation. In fact, this is how pulsars are 
observed here on Earth. The entire beam luminosity is 
collected as the signal sweeps across the telescope. Pur- 
suing this rough most conservative estimate for a pulsar 
at a distance d* = 100M, a signal 10~ 4 times fainter 
would be seen just behind the original pulse. The lumi- 
nosity is also transient, decaying in a timescale related to 
the instability of the orbits. (For a solar mass black hole 
r - 2GM Q /c 3 ~ 10" 5 - 1CT 6 sec.) Even if the beam 
pointed away from the Earth, as it swept over the black 
hole it would feed the stochastic layer and we would see a 
faint echo of the unseen pulsar from the diffuse light defo- 
cused off the black hole. The gravitational Doppler shift 
will also separate the frequency of the echo from that of 



the original pulse for nonstatic systems. As a last point, 
superradiance of scattered light from the ergosphere of a 
rotating black hole could also be significant. 

A more realistic calculation will be challenging as is re- 
flected by the imfamous full relativistic binary problem. 
What is clear is that as the companion gets closer, chaos 
will be more important and the signal will get brighter, 
but the lifetime of the binary will also be shorter. The 
last stages of inspiral may be characterized by the defo- 
cused echo in coordination with the gravity wave signal 
expected. 

Most black hole systems currently accessible to obser- 
vation involve the accretion of material from a luminous 
companion and the defousing effect would be completely 
obscured. Only in the most minimal binary pairs will 
the chaotic scattering lead to a visible glow around the 
black hole such as black hole/neutron star pairs and black 
hole/puslar pairs in particular. Since these are amoung 
the systems the future gravity wave experiments hope 
to discover, this could offer a valuable electromagnetic 
observational counterpart to any gravity wave detection. 
The gravity wave experiments hope to detect quite dis- 
tant coalescing binaries. It is unlikely that a pulsar would 
ever be visible at such large distances. Nonetheless, dur- 
ing the last stages of inspiral tidal stresses will undoubt- 
edly heat up the companion providing a brighter electro- 
magnetic signal to observe. The details of such a scenario 
are far from clear but the possibilities are worth investi- 
gating. 

While we have been promoting chaotic defocusing as 
a means to view the inner orbits around a black hole, 
there are other observable consequences of the chaotic 
flows. For instance gravity-waves are a natural and in- 
evitable source for the perturbations 0Jl^] and the im- 
plications for the direct detection of gravity-waves from 
the last stages of inspiral in a compact binary will cer- 
tainly be significant. More immediately observable, the 
disrupted motions of an accretion disk around black hole 
candidates could lead to an indirect detection of gravity- 
waves. Whether in radio waves or gravity waves, chaos 
may in fact make black holes bright. 
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[1] C.W. Misner, K.S. Thorne and J. A. Wheeler, Gravitation 
(W.H. Freeman and Company, New York, 1970). 

[2] C. Grebogi, E. Ott & J. A. Yorke, Phys. Rev. Lett. 50, 
935 (1983). 



5 



[3] E. Ott, Chaos in dynamical systems, (Cambridge Univer- 
sity Press, Cambridge, 1993). 

[4] S.D.Majumdar, Phys. Rev. 72 (1947) 390. 

[5] A. Papapetrou, Proc. R. Irish Acad. A51 (1947) 191. 

[6] G. Contopolous, Proc. R. Soc. A431 (1990) 183; G. Con- 
topolous, Proc. R. Soc. A435 (1990) 551. 

[7] C. P. Dettmann, N. E. Frankel and N. J. Cornish, Phys. 
Rev. D50, R618 (1994); Fractals, 3, 161 (1995). 

[8] N.J. Cornish and J.J. Levin, Phys. Rev. D 53 (1996) 
3022. 

[9] N.J. Cornish and J.J. Levin, Phys. Rev. Lett. 78 (1997) 
998; N.J. Cornish and J.J. Levin, Phys. Rev. D.55 (1997) 
7489. 

[10] N.J.Cornish and G.W. Gibbons, Class. Quantum Grav. 

14 (1997) 1865. 
[11] J.D.Barrow and J.Levin, Phys. Rev. Lett, 80 (1998) 656. 
[12] L. Bombelli and E. Calzetta, Class. Quantum. Grav. 9 

(1992) 2573. 

[13] R. Moeckel, Commun. Math. Phys. 150 415 (1992). 



G 



